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ABSTRACT 

We calculate theoretical population ratios of the ground fine-structure levels of some 
atoms/ions which typically exhibit UV lines in the spectra of QSO absorbers redward 
the Ly-ev forest: C°, C + , O , Si + and Fe + . The most reliable atomic data available 
is employed and a variety of excitation mechanisms considered: collisions with several 
particles in the medium, direct excitation by photons from the cosmic microwave 
background radiation (CMBR) and fluorescence induced by a UV field present. 

The theoretical population ratios are confronted with the corresponding column 
density ratios of C I and C n lines observed in damped Ly-a (DLA) and Lyman Limit 
(LL) systems collected in the recent literature to infer their physical conditions. 

The volumetric density of neutral hydrogen in DLA systems is constrained to be 
lower than tens of cm -3 (or a few cm -3 in the best cases) and the UV radiation field 
intensity must be lower than two orders of magnitude the radiation field of the Galaxy 
(one order of magnitude in the best cases) . Their characteristic sizes are higher than a 
few pc (tens of pc in the best cases) and lower limits for their total masses vary from 
10° to 10 5 solar masses. 

For the only LL system in our sample, the electronic density is constrained to be 
n e < 0.15 cm -3 . We suggest that the fine-structure lines may be used to discrimi- 
nate between the current accepted picture of the UV extragalatic background as the 
source of ionization in these systems against a local origin for the ionizing radiation 
as supported by some authors. 

We also investigate the validity of the temperature-redshift relation of the CMBR 
predicted by the standard model and study the case for alternative models. 

Key words: quasars: absorption lines - cosmic microwave background - atomic pro- 
cesses. 



1 INTRODUCTION 

Typically in the spectra of bright QSOs, several lines may be 
identified as beeing due to absorption of intervening material 
situated along the line of sight. These lines originate when 
the continuum radiation meets a common atom or ion in its 
ground state that partially absorbs it, leaving its imprint in 
the emiting QSO's spectrum. However, if some excitation 
mechanism is present in the absorbing region, then a small 
fraction of atoms or ions will also be found populated in their 
lowest-lying excited levels. Therefore, in addition to absorp- 
tion lines arising from the atom/ion's ground state, one may 
also expect to detect weaker lines arising from excited levels. 
It has been long pointed out that fine-structure absorp- 
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tion lines arising from the ground and lowest-lying excited 
energy levels of common atoms/ions may be used as an in- 
dicator of the physical conditions in the gas (Bahcall & Wolf 
1968; Smeding & Pottasch 1979). 

If we model the absorbing region as a single, homoge- 
neous cloud, then the ratio of the volumetric densities of 
atoms/ions populated in excited states n* to atoms/ions in 
the ground state n will match the corresponding column 
density ratios: 



For example, the column densities of C ions popu- 
lated in their ground 2 P°/ 2 and first excited 2 P-j/2 levels may 
be inferred from the equivalent widths of the corresponding 
2s 2 2p 2 P° /2 -> 2s2p 2 2 D° /2 and 2s 2 2p 2 P° /2 -> 2s2p 2 2 D° /2 
UV lines at 1334.5 A and 1335.7 A, respectively. 

The lefthand side of equation (1) in turn, may be evalu- 
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ated theoretically as a function of the physical conditions in 
the medium by solving the detailed equations of statistical 
equilibrium. It will in general depend on the intensities of 
several competing excitation mechanisms, such as sponta- 
neous decay, collisions with particles present in the medium 
or induced by radiation (the later either directly or by fluo- 
rescence) . 

The effectiveness of using column density ratios deduced 
from fine-structure lines to infer the basic parameters of a 
given excitation mechanism will depend on its relative im- 
portance to other processes contributing to the excitation 
of the fine structure levels. If collisions by a given particle 
dominate, one may expect to be able to infer its volumetric 
density; if fluorescence dominates, one is capable of measur- 
ing the intensity of the radiation field present, whereas if 
the dominating mechanism is direct excitation by photons 
of the Cosmic Microwave Background Radiation (CMBR) 
one could measure its temperature. This dominance of some 
process over another is determined not only by their relative 
intensities, but also by an interplay of atomic physics input 
parameters. 

The goal of this paper is to calculate theoretical popu- 
lation ratios of fine-structure levels of atoms/ions commonly 
found in QSO spectra, and to use them to estimate the phys- 
ical conditions in QSO absorbers. 

In section 2 we describe how the equations of statistical 
equilibrium were solved and details on the calculations for 
each selected atom/ion, namely: C°, C + , Si + , O and Fe + . 
In section 3 we gather recent column density ratios data 
taken from the literature and use the results obtained in 
the previous section to determine the physical conditions 
in Damped Lyman-a (DLA) and Lyman Limit (LL) QSO 
absorption line systems. For the DLA systems we also derive 
their characteristic sizes and masses. The main conclusions 
are sketched in section 4. 



2 ATOMIC PHYSICS 

In this section we calculate the population ratios of fine- 
structure levels for five atoms/ions of interest: C°, C + , 0°, 
Si + and Fe + . The first two already have their fine-structure 
lines observed with currently available high resolution spec- 
trographs (see section 4 below). The other atoms/ions are 
typically observed in QSO absorption spectra, since they 
have resonant lines longward the Ly-a line at 1216 A (so 
that they will not always fall into the Ly-a forest region of 
the spectrum). They also have their ground term split into 
fine-structure levels, and once the lines arising from excited 
levels are detected in future generations of more powerful 
telescopes, they may also provide useful physical condition 
information. 

Let us now briefly outline the basic procedures to cal- 
culate the fine-structure levels population ratios, and next 
discuss each particular atom/ion in greater detail. 

2.1 The statistical equilibrium equations 

In order to calculate the level populations of a given 
atom/ion, we make two basic assumptions : 

(i) The rates of processes involving ionization stages other 



than the atom/ion being considered (such as direct pho- 
toionization or recombination, charge exchange reactions, 
collisional ionization, etc.) are slow compared to bound- 
bound rates. 

(ii) All transitions considered are optically thin. 

In steady state regime the sum over all processes pop- 
ulating a given level i will be balanced by the sum over all 
processes depopulating it. Assuming that the two conditions 
listed above are met, this can be written: 

3 3 

where m is the volume density of atoms or ions in level 
i. We have defined the total rates taking the atom/ion from 
level i to level j as: 

Qij = Aij + BijUij + Yij + ^ nk< lij > ( 3 ) 

k 

where the coefficients Aij are transition probabilities, 
Bij are Einstein coefficients, Uij are the energy densities of 
the radiation field at the frequency of the transition i>ij, 
Fij are indirect excitation rates by fluorescence, n k are the 
volumetric densities of a given collision particle (usually k — 
e~ ,p + , H° , He , H2, ■ ■ ■ , depending whether the medium is 
primarily ionized or neutral) and qfj are the collision rates 
by some particle k. We have set Aij = for i < j and 
B u = u u = V li = = 0. 

Hereafter we shall abbreviate: 

Kij = BrjUij . (4) 

The indirect excitation rates are defined as (Silva & 
Viegas 2000): 

r = u ^mj + Kyi . . 

U 3 -^K^ {A ^ + R ^ , (5) 

pi v 

i.e., we have the situation in which the atom/ion - in 
one of its m lowest energy levels, i - is photoexcited to some 
higher energy level fj, and next decays - either spontaneously 
or by stimulated emission - back to some other level j among 
the lowest m. The sum extends over all possible upper levels. 

The fine-structure levels may also be directly populated 
by the CMBR. In that case one must add to the energy 
densities Uij the contribution from a black body radiation 
field redshited to a temperature (see, for instance, Kolb & 
Turner 1990): 

T = T (l + z) , (6) 

where T = 2.725 ± 0.001 K (la error) is the current 
value of the CMBR temperature as determined from the 
COBE FIRAS instrument (Mather et al. 1999; Smoot & 
Scott 2000). 

We caution, however, that this relation remains yet ob- 
servationally unprooven. In section 3.3 we review the cur- 
rently available pieces of evidence. 

Equation (2) is the system of statistical equilibrium 
equations that must be solved in order to compute the pop- 
ulation ratios of the fine-structure levels. If we model the 
atom/ion as beeing composed of m levels, then we must 
deal with a system of m — 1 equations. 

In order to numerically solve this system we have built 
a Fortran 90 code - popratio (Silva & Viegas 2000)- that 
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reads in the basic atomic physics parameters and automati- 
cally computes the rates for all the processes beeing consid- 
ered. The code is very flexible, allowing the user to account 
for an arbitrary number of levels and processes. 

The code, as well as the input files for the atoms/ions 
considered in this paper, are available upon request from 
one of the authors (AIS)"f. It may also be used in other as- 
tronomical applications, such as calculating intensity ratios 
of collisionally excited emission lines (such as coronal emis- 
sion lines) and computing cooling rates due to collisional 
excitation. 

Next, we describe the computations for each atom/ion 
considered in greater detail. As the population ratios of the 
fine-structure levels will be strongly dependent upon several 
atomic physics parameters, it is essential to search the liter- 
ature for the most up to date values. In this work, we give 
precedence to results obtained recently by two large inter- 
national collaborations: the Opacity Project (Seaton et al. 
1994) and the Iron Project (Hummer et al. 1993). 

For reasons of space, we illustrate the results obtained 
for the population ratios of the fine-structure levels under a 
limited range of physical conditions only. We urge the user 
to make use of the numerical code in order to get accurate 
predictions in his/her applications. 



2.2 The atom C° 

The ground state of the C° atom is comprised of the 2s 2 2p 2 
3 Po,i,2 triplet levels. The energies of the fine-structure ex- 
cited levels relatively to the ground state are 16.40 cm -1 
and 43.40 cm -1 . The transition probabilities are A w = 
7.932 10~ 8 s -1 , A 20 = 2.054 10~ 14 s" 1 and A 2 i = 
2.654 10~ 7 s" 1 . 

Our model atom includes the five lowest energy levels: 
2s 2 2p 2 3 Pf,,i, 2 , 2s 2 2p 2 and 2s 2 2p 2 ^g. The energies were 
taken from Moore (1970) and the transition probabilities 
from the Iron Project calculation of Galavi's, Mendoza and 
Zeippen (1997). 

The CMBR will be an important excitation mechanism 
for the first excited 3 P5 level, since it is so closely separated 
from the ground level. Assuming the temperature-redshift 
relation as given by equation (6), the CMBR spectrum will 
peak at the first-excited level frequency at a redshift z ~ 2. 
Table 1 gives the excitation rates of the C° fine-structure lev- 
els as a function of redshift, again assuming the temperature- 
redshift relation (6). 

The fine-structure transitions may also be induced by 
collisions. Fig. 1 shows the collision rates for the most im- 
portant collision particles. The rates for collisions by protons 
were taken from Roueff & Le Bourlot (1990), by neutral 
hydrogen from Launay & Roueff (1977a), by molecular hy- 
drogen from Schroder et al. (1991) and by neutral helium 
from Staemmler & Flower (1991). For the rates by collisions 
with electrons we have employed the analytic fits given by 
Johnson, Burke & Kingston (1987). We point out that the 
similar plot in Roueff & Le Bourlot's paper comparing col- 
lision rates by protons and electrons is incorrect, since an 
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Table 1. Excitation rates Kjj, of the C°, C+ and O fine- 
structure levels by the CMBR. We have assumed the temperature- 
redshift relation as predicted by the standard model (see text). 
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error has crept in their figure, and they compare excitation 
and de-excitation rates (Roueff, private comm.). 

We have also considered the effect of excitation of the 
upper 1 T>2 and 1 So levels by collisions with electrons. We 
took the analytic fits to the Maxwellian-averaged collision 
strengths 7 for transitions involving these levels given by 
Pequignot & Aldrovandi (1976), and transformed them from 
LS coupling to the fine-structure levels according to their 
statistical weights: 

7 (3pe j ^l D e ) = ^_ti 7 (3pc^l DC) (?) 
7 (3 P ^l S c) = ^i 7 (3pc^l SC) 

However, the inclusion of these levels can hardly influ- 
ence the population of the 3 P e fine-structure levels at tem- 
peratures prevailing in ionization regions where the atom C° 
is likely to be found. For example, even for temperatures as 
high as T — 10 4 K the population ratio of the 3 Pf level rel- 
atively to the ground state will increase by no more than 5 
percent (10 percent for the P§ level). The test calculations 
were done taking into account only collisions by electrons 
(and spontaneous decays) ; if this is not the main excitation 
mechanism, then the error will be significantly smaller. 

Excitation of the fine-structure levels by fluorescence 
was also investigated. We consider 108 allowed UV transi- 
tions involving the ground 3 P e levels and upper levels listed 
in the compilation of Verner, Verner & Ferland (1996), which 
is based on Opacity Project calculations. If we adopt the ra- 
diation field of the Galaxy (Gondhalekar, Phillips & Wilson 
1980), then the corresponding indirect excitation rates will 
be Toi = 3.5 10~ 10 s- 1 and r 02 = 2.8 10~ 10 s _1 . 

In fig. 2 we have plotted the population ratios of the C° 
fine-structure levels taking into account collisions by hydro- 
gen atoms (the main collision partner in ionization regions 
where the atom C° is found), the CMBR and fluorescence 
induced by the radiation field of the Galaxy. 

We compared our results with the previous calculations 
by Keenan (1989). He considered the effect of collisions by 
electrons and hydrogend atoms, as well as fluorescence in- 
duced by the radiation field of the Galaxy. Test calculations 
revealed good overall agreement with the values obtained by 
Keenan, with differences typically less than 15 percent. 

2.3 The ion C+ 

The ground state of the C + ion consists of the 2s 2 2p 2 P^ 3 

2 ' 2 

doublet levels. The energy of the fine-structure excited level 
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Figure 2. Population ratios of the C° fine-structure levels relatively to the ground state nj/riQ = n( 3 Pj)/n( 3 Pg) calculated under 
various physical conditions. The curve for the «,2/no population ratio for z = coincides with the curve taking only collisions into 
account. 



relatively to the ground state is 63.42 cm 1 , and the tran- 
sition probability is .43 1 = 2.291 1CT 6 s -1 . 

Our model ion includes the five lowest LS terms: 2s 2 2p 
2 P° and the 2s 2p 2 configurations 4 P C , 2 D C , 2 S C and 2 P C , 
making a total of ten levels when the fine-structure split- 
ting is accounted for. The energies were taken from Moore 
(1970) and the transition probabilities from the Iron Project 
calculation of Galavfs, Mendoza & Zeippen (1998). 

As the fine-structure levels of C + are more separated 
than the C° levels, the CMBR will play a significant role 
at higher redshifts only, as one can see from the excitation 
rates given in table 1 § . 

We take into account collisional excitation of the fine- 
structure levels with several particles. For the Maxwellian- 
averaged collision strengths for collisions by electrons we 
have adopted the calculation of Blum & Pradhan (1992). 
As their results differ by only 2 percent from the earlier cal- 
culation of Keenan et al. (1986), we have also included the 



§ Hereafter we shall assume as a working hypothesis the 
temperature-redshift relation predicted by the standard model. 



later's results at temperature values not covered by Blum & 
Pradhan's calculation as a means of broadening the avail- 
able temperature range. We took excitation rates by colli- 
sions with hydrogen atoms from Launay & Roueff (1977b), 
extrapolated to T > 1000 K by Keenan et al. (1986). Other 
collision particles taken into account are protons (Foster, 
Keenan & Reid 1997) and molecular hydrogen (Flower & 
Launay 1977). Fig. 3 compares the excitation rates with the 
various particles. 

We have complemented the work of Galavfs et al. with 
the allowed transitions listed in the compilation of Verner et 
al., making a total of 48 transitions involving the ground 2 P° 
levels and upper levels. The indirect excitation rate by the 
UV radiation field of the Galaxy could then be determined: 
Ti 3 = 9.3 10~ n s _1 . 

2 2 

In order to assess the relevance of the 2s 2p 2 configu- 
ration upper levels in the relative population of the ground 
2 P°i 3 levels, we have performed test calculations comparing 

2 ' 2 

our 10-level model ion with the 2-level ion. At high temper- 
atures the 2s 2p 2 configuration levels may be excited by 
collisions with hot electrons in the medium. However, the 
testcases have shown that this effect does not contribute 
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Figure 1. Excitation rates qjj/ = q( 3 Pj — > 3 P of the C° fine- 
structure levels by collisions with various particles. The points - 
taken from the literature cited in the text - are interpolated by 
cubic splines. 



Figure 3. Excitation rates qi 3 = <2( 2 -P? -* 2 -P|) of the C+ fine- 

22 2 2 

structure level by collisions with various particles. The points - 

taken from the literature cited in the text - are interpolated by 

cubic splines. 



significantly to the excitation of the 2 P° levels for tempera- 
tures T < 30000 K, where the discrepancies reach about 5 
percent. 

Therefore, for temperatures lower than 30000 K, only 
two levels can be taken into account. The population ratio 
of the excited fine-structure level relatively to the ground 
level is then expressed by: 



Ki 3 
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The collisional de-excitation rates may be computed 
from the principle of detailed balance: 



03 1 
^33 



-01 3 C 
2 33 



(9) 



with T expressed in K. 

In fig. 4 we have plotted the population ratio of the 
C + excited fine-structure level relatively to the ground state 
under various physical conditions. As the ion C + may coexist 
in both Hi and Hn regions, we sample two cases of interest: 
a neutral medium at T = 1000 K, and an ionized medium 
at T = 10000 K. In the later case, in addition to collisions 
by electrons, we also consider proton collisions and set n p — 
n e . However, at T — 10000 K their effect on the relative 
population ratio is only marginal (at the 5 percent level). 

Previous work on the population of the C + fine- 
structure levels taking into account fluorescence and colli- 
sions by electrons and hydrogen atoms was accoplished by 
Keenan et al. (1986). Test calculations showed that our re- 
sults seem to be in good agreement with their values, al- 
though it is not possible to make an accurate statement 
of the discrepancies, since they have published their results 
only in graphical form. 
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Figure 4. Population ratio of the C + fine-structure level rela- 
tively to the ground state nz/ni = n ( 2 P° )/n( 2 P° ) calculated 

under various physical conditions. The curves for z < 1 coincide 
with the curve taking only collisions into account. In the lower 
plot we have also taken proton collisions into account, n p = n e . 

2.4 The atom O 

The ground state of the O atom is comprised of the 2s 2 2p 4 
3 P 2 ,i,o triplet levels. The energies of the fine-structure 
excited levels relatively to the ground state are 158.265 
cm -1 and 226.977 cm~ . The transition probabilities are 
A12 = 8.865 1(T 5 s"\ A)2 = 1.275 1(T 10 s" 1 and A m = 
1.772 1(T 5 s _1 . 

Our model atom includes the five lowest energy levels: 
2s 2 2p 4 3 P 2 ,i, , 2s 2 2p 4 l T>% and 2s 2 2p 4 The energies were 
taken from Moore (1993) and the transition probabilities 
from the Iron Project calculation of Galavi's, Mendoza and 
Zeippen (1997). 

As the fine-structure levels of atomic oxygen are much 
more separated compared to atomic and singly ionized car- 
bon, the CMBR will not play a major role as one can see 
from the excitation rates for the first excited level given in 
table 1 (the excitation rates for the second excited level are 
even lower). 

The excited levels may be populated by collisions with 
particles present in the medium. Fig. 5 shows the collision 



rates for the fine-structure transitions induced by collisions 
with various particles. The rates for collisional excitation 
by electrons were taken from Bell, Berrington & Thomas 
(1998), by neutral hydrogen from Launay & Roueff (1977a) 
and by neutral helium from Monteiro & Flower (1987). For 
collisions with protons we have employed the analytic fits 
given by Pequignot (1990; 1996). 

For the sake of completeness, we have also considered 
collisional excitation of the upper X D| and 1 Sq levels. We 
have taken the Maxwellian-averaged collision strengths for 
transitions induced by electrons involving these levels from 
Berrington & Burke (1981). The rate for the Sp^D transi- 
tion induced by neutral hydrogen was taken from Federman 
& Shipsey (1983). The rates were transformed from LS cou- 
pling to the individual fine-structure levels according to eq. 
(7). 

After including 135 UV allowed transitions involving 
the ground 3 P C levels and upper levels from the work of 
Verner et al., we obtained the indirect excitation rates by 
the radiation field of the Galaxy: r 2 i = 3.9 1CP 11 s _1 and 
r 20 = 1.1 HT 11 s" 1 . 

The relative populations of the ground 3 PJ levels may 
be significantly affected by charge exchange reactions with 
hydrogen (Pequignot 1990; Pequignot 1996): 

H + + O ( 3 P° T ) -> H°( 2 Sl /2 ) + + ( 4 S° 3/2 ) (10) 
H°( 2 S5 /2 ) + 0+( 4 S° /2 ) -+ H+ + O'W) . 

Consideration of this process would require a knowledge 
of the ionization state of cloud, which lies beyond the scope 
of this paper. Therefore, in our analysis we consider only the 
case of a primarily neutral medium^. 

In fig. 6 we plot the population ratios of the ground 
O fine-structure levels under various physical conditions. 
We consider collisions by hydrogen and helium atoms, as- 
suming a helium abundance relative to hydrogen of 10 per- 
cent (by number). Collisions by helium atoms increases 
n( 3 P5)/n( 3 Pi) by only 5 percent and n( 3 PS)/n( 3 P^) by 10 
percent (reducing to zero close to LTE in the high density 
limit). The curves for n( 3 P5)/n( 3 P 2 ) corresponding to the 
inclusion of the effects of the CMBR at z — 5 and the UV 
field of the Galaxy are coincident because the relevant exci- 
tation rates are of the same order K^'° — ^21 ■ 

Our results are not directly comparable to the work 
of Pequignot (1990; 1996), who made assumptions on the 
ionization state of the gas. We point out, however, the im- 
portance of updating the electron excitation rates employed 
in their work - taken from Berrington (1988) - to the more 
recent calculations of Bell et al., since the later's results are 
substantially lower. 

2.5 The ion Si+ 

The ground state of the Si + ion consists of the 3s 2 3p 2 P^ 3 

2 ' 2 

doublet levels. The energy of the fine-structure excited level 
relatively to the ground state is 287.24 cm -1 , and the tran- 
sition probability is A3 1 = 2.17 10~ 4 s _1 . 

« This is the case for the DLA systems (section 3.1), where Ol 
lines arc commonly observed. 
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Figure 6. Population ratios of the O fine-structure levels relatively to the ground state nj/n,2 = n( 3 Pj)/n( 3 P|) calculated under 
various physical conditions. 



Our model ion includes the three lowest LS terms: 3s 2 3p 
2 P°, 3s 3p 2 4 P C and 3s 3p 2 2 D C , making a total of 7 levels 
when the fine-structure is accounted for. The energies were 
taken from Martin and Zalubas (1983). The transition prob- 
abilities for the 2 P^ — ► 2 P°i forbidden transition was taken 

2 2 

from Nussbaumer (1977), those for the 4 P e — > 2 P° intercom- 
bination transitions from Calamai, Smith & Bergeson (1993) 
and those for the 2 D e — > 2 P° allowed transitions from Nahar 
(1998). 

Since the fine-structure levels of Si + are too separated 
apart from eachother, the CMBR will not be an impor- 
tant excitation mechanism. Even for extremely high red- 
shifts z = 5, the excitation rate was found to be just 
K13 =4.7 1(T 15 s" 1 . 

2 2 

Collisional processes considered are collisions by elec- 
trons (Dufton & Kingston 1991), protons (Bely & Faucher 
1970) and hydrogen atoms (Roueff 1990). In fig. 7 we have 
plotted the excitation rates by collisions with these particles. 

Because the Maxwellian-averaged collision strength for 
the 2 P^ — 2 P*| transition induced by electrons varies by 

2 2 

no more than 6 percent in the calculated interval - 3.6 < 
logT < 4.6 - we also indicate in fig. 7 what might be ex- 



pected for the excitation rate down to T = 100 K if we 
assume a constant value for the collision strength. 

In order to account for fluorescence, we consider 39 UV 
allowed transitions from the work of Nahar (1998). The indi- 
rect excitation rate by the UV field of the Galaxy was found 
to be Pi 3 = 1.1 10~ 9 s" 1 . 

At sufficiently high temperatures the 4 P° and 2 D° up- 
per levels may be populated through collisions with hot elec- 
trons in the medium, and thereby influence the population 
of the P° ground levels. To assess the relevance of this ef- 
fect, we performed test calculations of the population ratios 
of the 2 P° fine-structure levels comparing the results ob- 
tained by the 2-level ion with those by the 7-level ion. Only 
collisions by electrons and spontaneous decays were consid- 
ered. The test calculations revealed that the upper levels 
are not important for T < 30000 K, when the discrepancies 
reach about only 6 percent. Therefore, as for C + , for tem- 
peratures lower than this, only two levels can be taken into 
account. The system of statistical equilibrium equations (2) 
then yields: 

n 3 Qis ri3+]r ,n k q k 13 

2_ _ 2 2 ry 2 2 ^ K 2 2 /-t-t\ 

in Q31 _ As i + V, n k q% 1 

2 2 2 2 2 L ~ lh ^l-J 
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Figure 5. Excitation rates qjji = q{ 3 Pj — * S Pji) of the O fine- 
structure levels by collisions with various particles. The points - 
taken from the literature cited in the text - are interpolated by 
cubic splines. 



Figure 7. Excitation rates qi 3 = q( P° — > Pf ) of the Si+ flne- 

22 2 2 

structure level by collisions with various particles. The points - 

taken from the literature cited in the text - are interpolated by 

cubic splines. The dotted line indicates an extrapolation of the 

excitation rate by electrons assuming a constant value for the 

corresponding Maxwcllian-avcragcd collision strength. 

Excitation and de-excitation collisional rates are related 

by: 



with T expressed in K. 

In fig. 8 we plot the population ratios of the fine- 
strutucture levels of Si + under various physical conditions. 
As Si + may be the prevailing ionization state in both H I 
and Hn regions, we sample two cases of interest: a neu- 
tral medium at T = 1000 K, and an ionized medium at 
T = 10000 K. 

Previous calculations of the population ratios of the 
ground fine-structure levels of Si + were performed by 
Keenan et al. (1985). Although test calculations under the 
same physical conditions considered by the authors ap- 
peared to reveal general agreement, it is difficult to quan- 
tify the discrepancies, since they published their results in 
graphical form only. We recommend the present calculations 
to the users, since they are based on more detailed and ac- 
curate atomic data. 

2.6 The ion Fe+ 

The ground state of the Fe + ion is comprised of the 

3d 6 4s 6 D e 9 7 5 3 1 sextet levels. Comparing to the other 

2 ' 2 ' 2 ' 2 ' 2 

atoms/ions previously studied, the ion Fe + has its fine- 
structure levels very separated apart from eachother and 
the transition probabilities are considerably higher. For ex- 
ample, the first excitated level is placed 384.790 cm -1 above 
the ground level, and the corresponding transition probabil- 
ity is ^79 = 2.13 10~ 3 s _1 . Both factors will contribute to 
make the population ratios of the fine-structure levels of the 
Fe + ion significantly low. 

Our model ion includes the four lowest LS terms: 3d 6 4s 
6 D C , 3d 7 4 F C , 3d 6 4s 4 D C and 3d 7 4 P°, making a total of 
sixteen levels when the fine-structure splitting is accounted 
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Figure 8. Population ratio of the Si + fine-structure level rela- 
tively to the ground state nz/ni = n ( 2 P° )/n( 2 P° ) calculated 

2 2 2 2 

under various physical conditions. In the dotted curves in the 
upper plot we have also added the contribution of fluorescence 
induced by the UV field of the Galaxy. 

for. The energies were taken from Corliss & Sugar (1982) and 
the transition probabilities from the Iron Project calculation 

of Quinet, Le Dourneuf & Zeippen (1996) H. 

Due to the high separation of the fine-structure levels, 
the CMBR will not be an important excitation mechanism. 
For example, at z — 5, the excitation rate to the first excited 
level is just K97 = 3.5 1CT 18 s" 1 . 

The only collisional process for which we could find de- 
tailed excitation rates calculated in the literature were col- 
lisions by electrons. Fig. 9 shows the excitation rates by 
collisions with electrons for the most important transitions 
within the 6 D e ground term. The corresponding Maxwellian- 



II We note that there is a small misprint for two transitions listed 
in table 5 of Quinet et al.'s paper. The authors did not add to 
the transition probabilitcs a magnetic dipole contribution, so that 
the values should actually read: a 4 F° 7 - a 4 P° 5 = 8.83 10~ 3 s _1 , 

a 4 F C 5 - a 4 P1 = 1.68 10 -3 s _1 , as it appears in their table 4. 
7 2 



1 3 10" 1 o 5 

- J'=7/2 




10 3 10" 10 5 

T(K) 

Figure 9. Excitation rates qg j; = q( 6 D% -> 6 D e ,) of the Fe+ 

2 2 

ground level to excited fine-structure levels by collisions with elec- 
trons. The corresponding Maxwcllian-avcragcd collision strengths 
were taken from the Iron Project calculation of Zhang & Pradhan 
(1995), and were interpolated by cubic splines. 



averaged collision strengths were taken from the Iron Project 
calculation of Zhang & Pradhan (1995). 

Nussbaumer & Storey (1980) estimated the excitation 
rates by collisions with protons to be less than 10 percent 
of the corresponding excitation rates by collisions with elec- 
trons for temperatures as high as T = 15000 K. However, 
as it is apparent from fig. 3 and fig. 7, the excitation rates 
for collisional processes involving positive ions and protons 
increase rapidly with temperature, so that one should be 
cautious when negleting collisions by protons at extremely 
high temperatures. 

We also include 212 allowed transitions involving the 
6 D e ground term levels and upper levels from the Iron 
Project calculation of Nahar (1995). The indirect excitation 
rates by fluorescence induced by the UV field of the Galaxy 
were found to be: Fg z = 7.0 10~ 10 s"\ Fg s = 1.3 10~ 10 s" 1 

2 2 2 2 

and To 3 = Tg i = (the later rates are zero because those 

2 2 2 2 

transitions are not eletric dipole allowed). 

Previous calculations of the fine-structure population 
ratios of Fe + levels were performed by Keenan et al. (1988). 
They included only the two lowest LS states in their model 
ion, 3d 6 4s 6 D° and 3d 7 4 F°, arguing that the next two LS 
states, 3d 6 4s 4 D e and 3d 7 4 P°, do not affect significantly the 
population of the 6 D e ground levels. However, test calcula- 
tions showed that considering the later LS terms increases 
the 6 D e ground level population ratios by as much as 17 
percent for T — 10000 K. The tests consisted of comparing 
the results of the 9-level and 16-level model ions taking only 
collisions by electrons (over various volume densities) and 
spontaneous decays into account. Therefore one should use 
the 16-level model ion for T > 10000 K. 

In order to assess the relevance at high temperatures of 
even higher-lying levels in the population ratios of the 6 D e 
ground levels, we have expanded our 16-level model ion to 
include the next two LS terms: 3d 7 2 G e and 3d 7 2 P C , thereby 
increasing the total number of levels to 20. 

Due to limitations of space, Quinet et al. list transition 
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Figure 10. Population ratio of the Fe + fine-structure levels rela- 
tively to the ground state nj/ng = n( 6 Dj)/n( 6 D c 9 ) as a function 

2 2 

of electronic density. 

probabilities just for the strongest (Aij > 10~ 3 s _1 ) transi- 
tions involving these levels. Hence, for the sake of complete- 
ness, we decided to complement their work with the weaker 
transitions from Garstang (1962). The Maxwellian-averaged 
collision strengths for these transitions were taken from the 
Iron Project calculation of Bautista & Pradhan (1996). 

Test calculations with the 20-level model ion revealed 
that the last two LS terms do not affect at all the popu- 
lation ratios of the 6 D° ground levels up to T = 20000 K 
(the highest temperature considered in Bautista & Prad- 
han 's calculation). 

Fig. 10 shows the population ratios of the Fe + fine- 
structure levels as a function of electronic density for T = 
10000 K. We may note a slight inversion in the population 
of levels 6 D 3 and 6 D1 at lower densities. 

2 2 

Comparison of our results with those from the previous 
calculation of Keenan et al. under the same physical con- 
ditions revealed that our values for the population ratios of 
the fine-structure ground levels are a factor of 2-4 larger. We 
believe this can be traced back to the Maxwellian-averaged 
collision strengths employed, since the values from Zhang & 
Pradhan are much higher than those obtained by Berring- 
ton et al. (1988), quoted by Keenan et al. Since the Iron 
Project calculation of Zhang & Pradhan delineates the reso- 
nance structure of the collision strengths in more detail, the 
results obtained in the calculation presented here should be 
more reliable. 



3 PHYSICAL CONDITIONS 

We now proceed to use our calculated atomic level popu- 
lation ratios to study the physical conditions in QSO ab- 
sorbers. 

Table 2 shows our sample of absorption line systems for 
which there are column density ratios of fine-structure lines 
reported in the recent literature. 

The sample includes DLA systems (logN(Hi) > 20.3), 
and only one LL system at z abs = 2.9034. 

We have not included any associated system, since their 



close proximity to the QSO could make them susceptible to 
the influence of the background radiation source, therefore 
requiring a case by case analysis that lies beyond the scope of 
this paper. The fine-structure lines are, however, a valuable 
tool to infer the physical conditions in such systems. In par- 
ticular, the knowledge of the ionization state of the systems 
coupled with the information on the volumentric density af- 
forded by the fine-structure lines allows one to place limits 
on the distance between the absorber and the QSO, giving a 
clue to infer whether they correspond to intervening clouds 
or to material ejected from the QSO (Turnshek, Weymann 
& Williams 1979; Morris et al. 1986; Tripp, Lu & Savage 
1996; Srianand & Petitjean 2000). 

So far, all the fine-structure lines observed belong to ei- 
ther C° or C + . Owing to its low ionization fraction (since its 
ionization potential is lower than that of hydrogen) , atomic 
carbon is very seldom detected. The three systems listed in 
table 2 correspond to all of the presently known C I systems, 
apart from the system observed towards the BL Lac object 
0215+015 (Blades et al. 1982; Blades et al. 1985). 

As we gathered observational data from the literature, 
we rejected any line falling within the Ly-a forest region 
of the spectrum. Prochaska (1999) observed the Cn* 1335 
fine-structure transition in a LL system at z a b s = 2.652 to- 
wards Q2231-00. However, since this transition falls within 
the Ly-a forest in this object and therefore may have been 
subject to significant contamination, his claimed value on 
the column density N(Cll*) should be regarded at most as 
an upper limit to the true value. For the same reason we dis- 
regarded the DLA system at z abs = 3.054 towards Q0000-26 
observed by Giardino & Favata (2000). Although the au- 
thors quoted their value for N(Cll*) as an upper limit, we 
argue that in principle significant contamination could also 
be taking place on the ground fine-structure line, thereby 
also affecting N(Cn) and driving the ratio N*/N in the op- 
posite sense. 

Unfortunately, the ground Cn 1334 line is often heavily 
saturated; to circunvent this problem there have been many 
alternative approaches to derive the N(Cn) column density 
by other indirect methods. Prochaska (1999) used the ratio 
of N(C n) /N(Fe n) in a velocity region where the ground C II 
line was not saturated to derive the corresponding value at 
the component where the C II* line was detected. Outram, 
Chaffee & Carswell (1999) assumed a carbon abundance rel- 
ative to iron [C/Fe]>-0.3 to obtain a tighter lower limit on 
the N(Cn) column density in a DLA system at z a b s = 2.62 
towards GB1759+75. In our sample we have included only 
direct measurements on the column densities. 

In sections 3.1-3.2 below, we will separately study the 
DLA and LL systems in our sample. Again, as a working 
hypothesis we shall assume the temperature-redshift relation 
as predicted by the standard model. The validity of this 
relation is discussed in section 3.3. 

3.1 DLA systems 

DLA systems have very high neutral hydrogen column densi- 
ties (log N(H i) > 20.3). This makes them effectively shielded 
from the ionizing radiation, causing their contents to be es- 
sentially neutral material (Viegas 1995). 

We use the fine-structure lines column density ratios 
observed in the DLA systems listed in table 2 to set upper 
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Table 2. Observational data on the column density ratios of fine-structure lines in QSO absorbers retrieved from the literature. 



# 


QSO 




^abs 


log N(H i) 


ion 


N*/N a 




^exc 


^cmbr 


reference 


1 


PKS 1756+23 


1.721 


1.6748 


> 20.3 


Ci 


< 1.2 10 _1 c 




< 7.4 


7.289 


Roth & Bauer (1999) 


2a 


Q1331+17 


2.084 d 


1.77638 


21.2 d 


Ci 


(3.1 ± 0.3) 10" 


-i 


10.4 ± 0.5 


7.566 


Songaila ct al. (1994b) 


2b 


Q1331+17 


2.084 d 


1.77654 


21.2 d 


Ci 


(1.3 ± 0.4) 10" 


-l 


7.4 ± 0.8 


7.566 


Songaila et al. (1994b) 


3 


Q0013-00 


2.0835 e 


1.9731 


20.7 ° 


Ci 


(4.0 ± 0.8) 10" 


-l 


11.7 ± 1.1 


8.102 


Ge, Bechtold & Black (1997) 


3 


Q0013-00 


2.0835 c 


1.9731 


20.7 ° 


Cn 


(7.0 ± 3.2) 10" 


■3 


16.2 ± 1.3 


8.102 


Ge, Bechtold & Black (1997) 


4 


Q0149+33 


2.43 


2.140 


20.5 f 


Cn 


< 9.6 10~ 3 




< 17.1 


8.557 


Prochaska & Wolfe (1999) 


5 


Q1946+76 


2.994 


2.8443 


20.27 


Cn 


< 2.1 10~ 2 




< 20.0 


10.476 


Lu ct al. (1996b) 


6 


Q0636+68 


3.174 8 


2.9034 


17.7 h 


Cn 


< 6.710~ 3 1 




< 16.0 


10.637 


Songaila et al. (1994a) 


7 


Q0347-38 


3.23 


3.025 


20.7 d 


Cn 


< 2.8 10~ 2 




< 21.4 


10.968 


Prochaska & Wolfe (1999) 


8 


Q2212-16 


3.992 


3.6617 


20.2 


Cn 


< 2.310- 2 




< 20.5 


12.703 


Lu ct al. (1996b) 


9 


Q2237-06 


4.559 


4.0803 


20.5 


Cn 


< 4.4 10~ 3 




< 14.9 


13.844 


Lu ct al. (1996b) 


10 


BRI 1202-07 


4.7 


4.3829 


20.6 


Cn 


< 1.2 10~ 2 




< 18.0 


14.668 


Lu ct al. (1996a) 



a Errors are la CL, while upper limits are 2a CL. 

b Assuming the temperature-rcdshift relation predicted by the standard model. 

c 2a upper limit on TV* obtained by private communication with the author. 

d Pettini ct al. (1994). 

° Ge & Bechtold (1997). 

f Wolfe et al. (1993). 

s Sargent, Steidel & Boksenberg (1989). 

h Derived from the optical depth of the LL discontinuity: tll = 3.5 (Sargent et al. 1989). 

1 Given the strong saturation of the ground fine-structure line, we adopt N > 1.510 14 cm~ 2 instead of the profile fitting value N = 
4.610 14 cm -2 preferred by Songaila ct al. (1994a). 



limits to their neutral hydrogen volume densities n H o and to 
the intensities of the UV radiation field present. Given the 
high neutral hydrogen column density, probably all of the 
hydrogen ionizing radiation will be absorbed, leaving very 
few photons with energies greater than 1 Ryd. The spectral 
shape of the UV radiation field will then be similar to the 
one found in our own galaxy, and we therefore assume the 
UV radiation field of Gondhalekar et al. multiplied by a 
constant factor /g- 

Table 3 shows the upper limits to n H o and fa for the 
DLA systems in our sample. They represent firm upper lim- 
its to the true values, because the single excitation mech- 
anism considered to obtain the upper limit - i.e., collisions 
by neutral hydrogen atoms to obtain n H o and fluorescence 
to obtain /g - may not be the dominating one and also be- 
cause for most systems the population ratios were just upper 
limits. 

Because the collisional excitation rate is temperature 
dependent, so will be the derived upper limits on n H o; wc 
assume two values of kinetic temperature characteristic of 
Hi regions: T = 100 K and T = 1000 K. 

If excitation by the CMBR is taken into account the 
upper limits become tighter (lower) , as indicated by the sec- 
ond figure next to each entry in table 3 (if only one value 
appears, it remains unchanged to the last significant digit 
displayed). Accounting for the CMBR affects considerably 
the results for the Ci systems (objects 1 and 2), since it 
is an important excitation mechanism for C° as mentioned 
earlier in section 2.2. Note the striking difference between 
both values for object 1, which has an excitation tempera- 
ture very close to the predicted CMBR temperature. We can 
not consider the CMBR for object 2b, since the excitation 
temperature is slightly lower than the CMBR temperature. 
As for the remaining C n systems the result is changed signif- 
icantly only for the z > 4 regions (objects 9 and 10), when 
the CMBR starts to play a significant role at the neutral 
hydrogen densities involved (cf. fig. 4, top). 



Collisions with molecular hydrogen are not likely to be 
relevant in our analysis, since the molecular fraction usu- 
ally seen in DLA systems is exceedingly small: /(Eh) = 
2N(H 2 )/N(H) < 210~ 4 , reaching as low as /(H 2 ) = 4 HT S 
in the z abs = 3.3901 DLA system towards Q0000-26 (Lev- 
shakov et al. 2000). Two exceptions are object 3 in our sam- 
ple (Ge & Bechtold 1997) and the z abs = 2.34 DLA system 
towards Q1232+08 (Ge, Bechtold & Kulkarni 2000), with 
/(Eh) = 0.22 and 0.07, respectively. In any case that would 
imply riH 2 g" 2 << n H ogg , as typically qf 2 < (cf. figs. 
1 and 3). 

From table 3 we see that the ratio of fine-structure lines 
observed in DLA systems constrain their neutral hydrogen 
densities to be lower than tens of cm~ 3 (or a few cm~ 3 in 
the best cases), and the UV radiation field to be lower than 
two orders of magnitude the radiation field present in our 
galaxy (or one order of magnitude in the best cases). 

Naturally, we could also have placed upper limits to the 
electron density n e . The upper limits on n e derived from C n 
lines would be about two orders of magnitude lower than the 
corresponding upper limits on n H o listed in table 3, i.e., in 
the approximate inverse ratio of the corresponding collision 

rates <7i 3/513 ~ 10 2 (fig. 3). For C° this ratio is no more 

22 22 

than 10 in the relevant temperature region (fig. 1). As DLA 
systems are comprised of mostly neutral material, the free 
electrons will come mainly from neutral atoms which have 
an ionization potential lower than that of hydrogen, such 
as C° , an whose (solar) elemental abundance relative to hy- 
drogen is of order 10~ 4 . Therefore, we would expect before- 
hand n e ~ 10~ 4 n H o, and the fine-structure lines would not 
provide a meaninful constraint on the electron density. The 
electron/neutral hydrogen density ratio may be even lower 
if we consider that DLA systems may exhibit abundances as 
low as two orders of magnitude below solar (e.g. Pettini et 
al. 1994). 

We can also estimate the characteristic sizes and to- 
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Table 3. Physical conditions in DLA systems. The quoted values for /q and n H o are upper limits, 
whereas those for I and M are lower limits. The CL is 2<r. The second figure next to each entry 
corresponds to the inclusion of the CMBR (if only one value appears, it is not altered to the last digit 
displayed) . 









T = 100 K 






T = 1000 K 




# 


/g 


n H o [cm 3 ] 


I [pc] 


M [M ] 


n H o [cm 3 ] 


I [pc] 


M [M Q ] 


1 


16 / 0.79 


16 / 0.79 


4.0 / 81 


26 / 10600 


6.9 / 0.34 


9.4 / 192 


142 / 58800 


2a 


42 / 25 


43 / 26 


11 / 19 


1540 / 4290 


18 / 11 


27 / 46 


8910 / 24900 


2b 


16/- 


16/- 


30 / - 


10700 / - 


6.9 / - 


70 / - 


59600 / - 


4 


236 / 235 


35 / 34 


2.8 


19 


12 


7.8 


146 / 148 


5 


518 / 510 


77 / 75 


0.79 / 0.80 


0.93 / 0.96 


27 


2.2 


7.3 / 7.5 


7 


703 / 691 


105 / 103 


1.6 


10 


37 / 36 


4.4 


77 / 79 


8 


582 / 544 


86 / 81 


0.60 / 0.64 


0.45 / 0.52 


31 / 29 


1.7 / 1.8 


3.6 / 4.1 


9 


107 / 39 


16 / 5.7 


6.2 / 17 


93 / 692 


5.7 / 2.1 


17 / 47 


711 / 5320 


10 


307 / 209 


45 / 31 


2.9 / 4.2 


26 / 57 


16 / 11 


8.0 / 12 


205 / 445 



tal masses of the absorbing clouds responsible for the DLA 
systems; these will be given by: 



I 

M 



N(H) 
n H 

m p nnl =m p N(H)l 



(13) 



where m p is the proton's mass and N(H) and riH are the 
total hydrogen column and volume densities, respectively. 

For DLA systems we can make the replacement N(H) = 
N(Hl) and «h — n u o ** ; hence we can use our upper limits 
on n H o to set lower limits to the characteristic sizes and 
total masses of the intervening clouds. We see from table 3 
that our constraints imply characteristic sizes larger than a 
few pc (tens of pc in the best cases) and lower limits for the 
total masses that vary from 10° to 10 5 solar masses. 

In deriving the cloud sizes and masses above, we have 
implicitly assumed that most of the hydrogen column den- 
sity is in the same component where the fine-structure lines 
could be measured. Altough it would be much more difficult 
to detect the excited fine-structure line in the velocity com- 
ponent with the lowest associated hydrogen column density, 
we can not rule out the possibility that this is compensated 
by a higher metalicity and intensity of local excitation mech- 
anisms. 

We now focus our attention to object 3 in our sam- 
ple, which exhibits both Cl and Cn fine-structure lines. In 
fig. 11 we derive the neutral hydrogen volume density as 
a function of the intensity of the UV radiation field based 
upon the column density ratios of C I and C n lines (and 
for the two values of kinetic temperature considered before; 
the CMBR is included). If we assume that C° and C + are 
located within the same ionization region in the cloud, then 
the physical conditions will be described by the intersection 
of both curves. For T = 100 K we have n H o = 24 cm~ 3 and 
fa — 11, whereas for T = 1000 K we have n H o = 14 cm~ 3 
and /g = 8. Therefore, regardless of the kinetic temperature 
adopted, the UV field present must be one order of magni- 
tude more intense than in our galaxy. This contrasts with 
object 1, where the observed C I lines constrain the UV field 



** Some authors use the C II fine-structure lines to constrain n e 
and set njj — n e ; from the discussion in the preceding paragraph 
we note that this underestimates «h by two orders of magnitude. 
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Figure 11. Physical conditions in object 3 in our sample (the 
CMBR is included). 



to be lower than in our galaxy. In the detailed photoioniza- 
tion model constructed by Ge et al. the physical conditions 
prevailing in most regions of the cloud are: T = 100 K, 
n H0 = 21.0 ± 9.6 cm" 3 , f G = 17.0 and n e = 5.010" 4 n H0 . 
Assuming their value for n H o we have I — 7.7 ± 3.6 pc and 



M = 240 ± 160 M, 



©• 



3.2 LL systems 

The LL systems differ considerably from the DLA systems 
studied before for beeing significantly ionized. The source of 
the ionizing radiation in these systems is usually assumed 
to be the UV extragalactic background, as the integrated 
radiation field of all QSOs attenuated by the intergalactic 
medium (Haardt & Madau 1996). 

Some authors claim for a local origin to the source of 
ionization (Viegas & Friaga 1995). In particular, the puz- 
zling observations of column densities of C, N and O ions in 
various ionization states as well as of He I in LL systems to- 
wards Q1700+64 (Vogel & Reimers 1993; Reimers & Vogel 
1993) cannot be simultaneously explained by photoioniza- 
tion models based on the UV background as the source of 
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Table 4. Indirect excitation rates of C + fine-structure levels by 
the radiation fields predicted by hot halo models. Each model, 
taken from Vicgas & Friaga, corresponds to a given age and dis- 
tance from the center of the forming galaxy. 



t (Gyr) r (kpc) 



Ti 3 

2 2 



(s-i) 



0.206 
0.206 
0.3644 
0.3644 



30 
100 
30 
100 



6.8 10 

6.2 10 

1.3 10 
1.2 10 



-10 

-ii 
-ii 

-12 



ionization. These systems are successfully interpreted by the 
hot halo model (Viegas & Friaga 1995); in this model the LL 
systems are identified as cold condensations embedded in a 
hot halo formed during the early stages of galaxy evolution, 
which acts as the source of ionization. 

Photoionization models typically constrain the ioniza- 
tion parameter U = ^a. (<j> H ; s the total photon flux of 
hydrogen ionizing radiation); the fine-structure lines might 
be used to independently constrain the volumetric density 
and test the hypothesis of a given radiation field as beeing 
the source of ionization. If the density determined from the 
fine-structure lines, together with the ionization parameter 
determined from ionization models based on the UV back- 
ground imply a higher intensity for the ionizing radiation 
field than expected, that would be a strong evidence for a 
local origin to the true ionization source. 

For the only LL system in our sample (object 6 in table 
2), we derive a 2a upper limit to the electronic density of 
n e < 0.15 cm -3 , assuming a kinetic temperature T = 10 4 K 
characteristic of photoionized regions. We have included the 
minor contribution from the CMBR and collisions by pro- 
tons (assuming a fully ionized medium n p = n e ), although 
they affect the result only at the 10 percent level. Fluores- 
cence plays a negligible role. In table 4 we show the indirect 
excitation rates of C + fine structure levels for the hot halo 
models considered by Viegas & Friaga; in any case we have 
r i 3 << n e g e i 3 . The indirect excitation rate induced by the 

2 2 2 2 

UV background turns out to be even lower; we have adopted 
the revised calculation of Madau, Haardt & Rees (1999) to 
obtain an indirect excitation rate at the observed redshift of 



3.3 The CMBR temperature-redshift relation 

The CMBR constitutes one of the cornerstones of the hot 
Big Bang model, which makes three basic quantitative pre- 
dictions on its properties: 

(i) it is isotropic and homogeneous; 

(ii) it has a black-body spectrum, 

(iii) it cools as the universe expands according to the re- 
lation T = T (l + z). 

Over the past decade, the advent of the COBE satellite 
has allowed the confirmation of the isotropy (Smoot et al. 
1992) and black-body spectral shape (Mather et al. 1994) 
to unprecedented precision, giving a present day tempera- 
ture of To = 2.725 ± 0.001 K (la error) as determined from 
the FIRAS instrument (Mather et al. 1999; Smoot & Scott 
2000). 



Any direct means of measuring the CMBR temperature 
can just provide us with the value of its current tempera- 
ture, forcing us to resort to other indirect methods to test 
the temperature-redshift relation predicted by the standard 
model. The best alternative is to use atomic and molecu- 
lar transitions seen in the spectra of QSO absorbers (Meyer 
1994). It is worth noting that the observation of CN ab- 
sorption lines from diffuse interstellar clouds towards bright 
stars in the Galaxy yielded T = 2.729^° q 3 j K, in excellent 
agreement with the COBE FIRAS result (Roth 1992; Roth, 
Meyer & Hawkins 1993; Roth & Meyer 1995). 

Unfortunately, molecular transitions are not commonly 
seen in the spectra of QSO absorbers. Apart from H2, so far 
molecules have been identified in just four absorption sys- 
tems (Wiklind & Combes 1994; Wiklind & Combes 1995; 
Wiklind & Combes 1996a; Wiklind & Combes 1996b). Sur- 
prisingly, in one of them (Wiklind & Combes 1996b) the ro- 
tational transitions from several molecules indicated an exci- 
tation temperature T oxc = 4±2 K (3<r error), lower than the 
expected CMBR temperature T — 5.14 K predicted at the 
observed redshift. The low excitation temperature in this 
object is, nevertheless, due to the effect of a microlensing 
event (Combes 2000, private comm.). Molecular absorption 
systems are often gravitational lenses, since the impact pa- 
rameter to the foreground galaxy must be close to zero in 
order to allow the detection of molecules. Hence, we believe 
that atomic lines are better suited to study the temperature 
of the CMBR at high redshifts. 

We can use the population ratios of the fine-structure 
levels for the absorption systems collected in table 2 to con- 
strain the temperature of the CMBR at their redshifts. For 
each observed ion the excitation temperature will be given 
by 



23.60 



ln(3/(n*/n)) 
91.25 



for C° 



for C 4 



(14) 



ln(2/(n*/n)) 
with temperatures given in K. 

The excitation temperatures so obtained represent firm 
upper limits to the temperature of the CMBR, because lo- 
cal excitation mechanisms may also contribute significantly 
to populate the excited levels. In fig. 12 we plot the excita- 
tion temperatures along with the expected temperature of 
the CMBR according to the standard model prediction. For 
most systems, either the signal to noise ratio of the spec- 
trum was not high enough to detect the excited fine struc- 
ture line, or the ground C n line was strongly saturated. 
Therefore for these systems the excitation temperature it- 
self is also an upper limit, and this is indicated in fig. 12 
by a downward arrow. The point labelled "molecules" cor- 
responds to the puzzling observation of Wiklind & Combes 
(1996b) discussed above. 

Phillips (1994) supports a closed steady-state model 
that predicts considerably lower temperatures to the CMBR 
compared to the standard model; e.g. for the z = 2.9 Cn 
system observed by Songaila et al. (1994a): T obs = 6.55 K. 

Alternative models in which photon creation takes place 
as the Universe expands predict a more general temperature- 
redshift relation (Lima, Silva & Viegas 2000): 

T = T„(l + z) 1 ~' 3 , (15) 
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Figure 12. Excitation temperatures derived from fine-structure 
absorption lines. The solid line is the temperature of the CMBR 
according to the tcmpcraturc-rcdshift relation given by the stan- 
dard model, while alternative models with photon creation predict 
a lower temperature (dotted lines). Error bars are la CL, whereas 
upper limits are 2a CL. 



where j3 is a parameter to be adjusted from the observa- 
tions, within the range < j3 < 1. Equation (15) therefore 
gives temperatures lower than predicted by the standard 
model (fig. 12). 

It has been stated that any scenario that does not pre- 
serve the number of photons would introduce large distor- 
tions in the black-body spectrum of the CMBR (Steigman 
1978). However, it was shown that for the class of models 
that follow the temperature law (15) the Planckian spectral 
shape is not destroyed as the Universe evolves (Lima 1996; 
Lima 1997). In these models the total entropy of Universe 
increases with time, but the entropy per particle remains 
constant. 

Big Bang nucleosynthesis arguments, however, severely 
limit the value of the free parameter to f3 < 0.13 (Birkel & 
Sarkar 1997). 

Inspection of fig. 12 reveals that current measurements 
do not require any extra ingredients to the standard model, 
since the totality of the points lie above the linear temper- 
ature law. However, a conclusive statement could only be 
made after correcting for local excitation mechanisms, in 
order to convert the excitation temperature upper limits to 
the actual temperature of the CMBR. Altough some points 
in fig. 12 appear to be dangerously close to the standard 
model prediction, it could be that additional local excita- 
tion mechanisms are negligible compared to excitation by 
the CMBR in these systems. In that case, the excitation 
temperature would provide a direct measure of the CMBR 
temperature. 

For object 3 in our sample, Ge et al. (1997) con- 
structed a detailed photoionization model to account 
for the local excitation mechanisms. They obtained 
T = 7.9 ± 1.0 K, whereas the standard model prediction is 
T — 8.102 ± 0.003 K (fig. 12). 



If the temperature law given by the standard model 
turns out to be incorrect, it would pose a serious source 
of difficulty, since not even the presence of a cosmological 
constant would alter the predicted temperature-redshift re- 
lation (Lima, Silva & Viegas 2000). On the other hand, if 
it is confirmed by the observations, that would add another 
success to its list of triumphs with a bonus: since each ab- 
sorbing region is located at a different site of the Universe, 
we could also assess its homogeneity. 



4 CONCLUSIONS 

We have presented new theoretical calculations of popula- 
tion ratios of the ground fine-structure levels of C°, C + , O , 
Si + and Fe + . The literature was searched for the most recent 
and reliable atomic data available to date. Various possible 
excitation mechanisms are taken into account. We encour- 
age the user to make use of the available Fortran code to 
obtain accurate predictions in his/her applications, rather 
then rely on approximate practical formulas that are valid 
only for a limited range of physical conditions. 

We have retrieved from the literature observational data 
on the column density ratios derived from the fine-structure 
lines, and confronted them with our theoretical calculations 
to infer the physical conditions prevailing in DLA and LL 
systems. Currently only C I and C n fine-structure lines have 
been observed; future detection of lines originating from less 
excited atoms/ions such as O , Si + and Fe + (which also have 
resonant lines redward the Ly-a forest) might aid to better 
constrain the physical conditions. 

For the DLA systems, the neutral hydrogen volumetric 
density is lower than tens of cm" 3 (a few cm" 3 in the best 
cases) and the UV radiation field is less intense than two 
orders of magnitude the UV field of the Galaxy (one order 
of magnitude in the best cases). Their characteristic sizes 
are higher than a few pc (tens of pc in the best cases) and 
lower limits for their total masses vary from 10° to 10 5 solar 
masses. 

For the only LL system in our sample, we derived n e < 
0.15 cm" 3 . As more observations become available, it may 
be possible to use the information contained in the fine- 
structure lines to help determine the nature of the source of 
ionization of these systems. 

The fine-structure lines in QSO absorbers also provide 
a method to test the temperature-redshift relation for the 
CMBR predicted by the standard model. Current observa- 
tions do not contradict the linear temperature law, although 
a conclusive statement could only be made after account- 
ing for local excitation mechanisms in order to correct the 
excitation temperatures to the actual temperature of the 
CMBR. That would require a knowledge of the ionization 
state of the cloud, after appropriate modeling by a photoion- 
ization code. 

A substantial improvement from the theoretical stand- 
point could be achieved by analysing together the ionization 
state of the cloud and the excitation of the fine-structure lev- 
els, by coupling our code - popratio - to a larger photoion- 
ization code. Presently, all studies based on the excitation of 
the fine-structure levels were carried out separately from the 
photoionization modeling, considering average physical con- 
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ditions throughout the entire cloud (e.g. Ge et al., Giardino 
& Favata). 

By carrying out both analyses simultaneously, we could 
further refine our models and eliminate the need for all 
assumptions we have maden in our analyses of the fine- 
structure lines. The intensities of the excitation mechanisms 
could vary across the cloud, and it would no longer be neces- 
sary to stick to the simplistic case of a perfect homogeneous 
cloud, that we implicitly assumed when we wrote down eq. 
(1). Moreover, we could generalize our statistical equilibrium 
equations (2) to include terms which require the knowledge 
of the ionization state of the cloud, such as recombination 
and charge exchange. It would also be possible to account 
for optical depth effects, allowing us to drop the assumption 
that all transitions considered are optically thin. 

On the observational side, there is also clearly a need 
for better measurements, since for the great majority of the 
systems only upper limits to the column density ratios are 
available. As the future generations of more powerful tele- 
scopes equipped with high resolution spectrographs continue 
to push the detection limits to even weaker lines, more in- 
formation could be available by observing atoms/ions other 
than C° and C+. 
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